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Abstract 

We study in this work the two-body hadronic charmed meson decays, including both the PP and 
VP modes. The latest experimental data are first analyzed in the diagrammatic approach. The 

^ ; 

magnitudes and strong phases of the flavor amplitudes are extracted from the Cabibbo-favored 
, (CF) decay modes using x 2 minimization. The best-fitted values are then used to predict the 

branching fractions of the singly-Cabibbo-suppressed (SCS) and doubly-Cabibbo-suppressed de- 
cay modes in the flavor SU(3) symmetry limit. We observe significant SU(3) breaking effects in 
some of SCS channels. In the case of VP modes, we point out that the Ap and Ay amplitudes 
cannot be completely determined based on currently available data. We conjecture that the quoted 
experimental results for both Df — > K°K* + and Df — > p + r]' are overestimated. We compare the 
sizes of color-allowed and color-suppressed tree amplitudes extracted from the diagrammatical ap- 
proach with the effective parameters a\ and a<i defined in the factorization approach. The ratio 
|a2/ai| is more or less universal among the D — > Kit, K*tt and Kp modes. This feature allows 
us to discriminate between different solutions of topological amplitudes. For the long-standing 
puzzle about the ratio T{D° K + K~)/T(D° — > 7r + 7r ), we argue that, in addition to the SU(3) 
breaking effect in the spectator amplitudes, the long-distance resonant contribution through the 
nearby resonance /o(1710) can naturally explain why D° decays more copiously to K + K~ than 
7r + 7r~ through the IT-exchange topology. This has to do with the dominance of the scalar glueball 
content of /o(1710) and the chiral-suppression effect in the decay of a scalar glueball into two 
pseudoscalar mesons. The same FSI also explains the occurrence of D° — > K°K° and its vanishing 
amplitude when SU(3) flavor symmetry is exact. Owing to the G-parity selection rule, Df — > tt + uj 
does not receive contributions from the short-distance IT-annihilation and resonant FSIs, but it 
can proceed through the weak decays D+ — > p + rj^ followed by the final-state rescattering of p + ij^ 
into tt + uj through quark exchange. 
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I. INTRODUCTION 



The hadronic decays of charmed mesons and related physics have been studied extensively in the 
past 35 years since the first discovery of charmed meson states in 1974. The experimental progress 
is still impressive. For example, many hadronic multibody decays of D°, D + and mesons have 
been studied using the technique of Dalitz plot analysis in which the resonant structure is probed. 
The Dalitz analysis offers the opportunity to understand the light scalar meson spectroscopy and 
the underlying structure of light scalar mesons a, K,ao(980) and /o(980). The study of charm 
physics also opens a new avenue to the search of new physics. For example, the recent observation 
of D° — D° mixing will enable us to explore new physics via flavor-changing neutral currents in 
the up-type quark sector. The standard model predicts almost null CP asymmetries in the charm 
system. Hence, any observation of CP violation in hadronic charm decays will signal physics beyond 
the standard model. 

With the advent of heavy quark effective theory, it is known that physics of nonleptonic B 
decays can be formulated in a QCD-inspired approach such as QCD factorization, pQCD and 
soft-collinear effective theory. However, until today a theoretical description of the underlying 
mechanism for exclusive hadronic D decays based on QCD is still not yet available. This has to do 
with the mass of the charm quark, of order 1.5 GeV. It is not heavy enough to allow for a sensible 
heavy quark expansion and it is not light enough for the application of chiral perturbation theory. 
It is known that the naive factorization approach fails to describe color-suppressed (class-II) decay 
modes. Empirically, it was learned in the 1980s that if the Fierz-transformed terms characterized 
yy 1/N C are dropped, the discrepancy between theory and experiment will be greatly improved 
l| . This leads to the so-called large- N c approach for hadronic D decays 0] . Theoretically, explicit 
calculations based on the QCD sum-rule analysis Q] indicate that the Fierz terms are indeed largely 
compensated by the nonfactorizable corrections. However, the 1/N C ansatz is not applicable to B 
decays. 

In spite of the difficulties in formulating a theoretical framework for describing nonleptonic 
decays of charmed mesons, a model-independent analysis of charm decays based on the diagram- 
matic approach is possible. In this approach, the flavor-flow diagrams are classified according to 
the topologies of weak interactions with all strong interaction effects included. Based on flavor 
SU(3) symmetry, this model-independent analysis enables us to extract the topological amplitudes 
and sense the relative importance of different underlying decay mechanisms. When enough mea- 
surements are made with sufficient accuracy, we can extract the diagrammatic amplitudes from 
experiment and compare to theoretical estimates, especially checking whether there are any sig- 
nificant final-state interactions or whether the weak annihilation diagrams can be ignored as often 
asserted in the literature. 

A salient feature of the diagrammatic approach analysis of charm decays is that, unlike the 
B decays, the magnitude and relative strong phase of each individual topological diagram can 
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be extracted from the data without making further assumptions apart from flavor SU(3) sym- 
metry. Explicitly, the color-allowed amplitude T, the color-suppressed amplitude C, TV-exchange 
amplitude E and VF-annihilation amplitude A in Cabibbo-allowed D — > PP decays can be cleanly 
extracted from the data. Recall that the ratio \C/T\ in charmless B decays cannot be unambigu- 
ously determined due to the contamination of electroweak penguin or weak annihilation diagrams 
and that weak annihilation topologies in the B sector are usually neglected in the diagrammatic 
approach analysis. 

The topological diagram analysis of D decays evidently shows the importance of weak annihi- 
lation (VF-exchange or TF-annihilation) contributions. Since the short-distance weak annihilation 
amplitude is subject to helicity suppression, this implies that long-distance weak annihilation in- 
duced by final-state interactions should play an essential role. 

In this work we shall study the two-body hadronic charmed meson decays, including both the 
PP and VP modes. We first analyze the latest experimental data in the diagrammatic approach. 
For D — > PP decays, there are new CLEO measurements with better precision Q]. There are two 
different ways of extracting the topological amplitudes for D —> VP decays. We shall present both 
results and make a comparison between them. At present, there exist many possible solutions for 
the flavor amplitudes. We will study the effective parameters a\ and 02 defined in the factorization 
approach and use them as a means to discriminate among different solutions. Another goal of this 
work is to see if we can understand the SU(3) breaking effects occurring in the SCS modes, in 
the hope of shedding more light on or even resolving the long-standing puzzle with the ratio of 
D° — > K + K~ to D° — > 7r + 7r~ branching fractions. 

The layout of the present paper is as follows. After a brief review of the diagrammatic approach, 
we extract magnitudes and strong phases of the flavor amplitudes from CF D — > PP and D — > VP 
decays in Sec. II. From the extracted invariant amplitudes, we determine the effective parameters 
a\ and 02 defined in the factorization approach in Sec. Ill and proceed to discuss SU(3) breaking 
effects in singly-Cabibbo-suppressed modes. Using the branching fraction ratio of D° — > K + K~ 
to D° —J- 7r + 7r~ as an example, we illustrate how to explain the large SU(3) violation in the ratio 
by considering the PF-exchange topology induced by the nearby resonance. Sec. IV comes to our 
conclusions. 

For nice and comprehensive reviews of exclusive hadronic charmed meson decays, the reader is 
referred to Refs. [3, M. 
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II. DIAGRAMMATIC APPROACH 



A. Preamble 

It has been established sometime ago that a least model-dependent analysis of heavy meson 
decays can be carried out in the so-called quark-diagram approach . In this diagrammatic scenario, 



all two- body non 
quark diagrams 



eptonic weak decays of heavy mesons can be expressed in terms of six distinct 
7H9(]: 2 T, the color-allowed external IF-emission tree diagram; C, the color- 
suppressed internal VF-emission diagram; E, the IF-exchange diagram; A, the ^-annihilation 
diagram; P, the horizontal IF-loop diagram; and V , the vertical VF-loop diagram. (The one-gluon 
exchange approximation of the P graph is the so-called "penguin diagram".) It should be stressed 
that these diagrams are classified according to the topologies of weak interactions with all strong 
interaction effects encoded, and hence they are not Feynman graphs. All quark graphs used in 
this approach are topological and meant to have all the strong interactions included, i.e., gluon 
lines are included implicitly in all possible ways. Therefore, analyses of topological graphs can 
provide information on final-state interactions (FSIs). Various topological amplitudes in two-body 



hadronic D decays have been extracted from the data in Ill4l8j] after making some reasonable 
approximations, e.g., flavor SU(3) symmetry. 

For charm decays involving an SU(3) singlet in the final product, e.g., D° — > ~K°(j), ~K°lo, K°rjo, 
there exist additional hairpin diagrams in which a quark-antiquark pair is created from vacuum 
to form a color- and flavor-singlet final-state meson [l{]; 13 • There are four different types of 
disconnected hairpin diagrams: Eh, Ah, Ph, Vh corresponding to the topological graphs E, A, P, V 



(for details, see [13]). The amplitudes Eh and Ah are denoted by SE and SA, respectively, in 



Ref. [121 ]. They appear in final states with rj^ for the PP modes and with (f>, to, and/or rf'> for the 



PV modes. A test of their importance has been proposed in [12|. As they are still not strongly 
called for phenomenologically (see discussions in Section III Bp , we will ignore these diagrams in 
the current analysis. For hadronic charm decays, contributions from the diagrams P and V are 
negligible due to the cancellation of the CKM matrix elements V* d V u d and V* S V US . 

For final states involving an r/ or rj' , it is more convenient to consider the flavor mixing of r\ q 
and rj s defined by 

Vq = ^( M ^ + = SS, (1) 

in analogy with the wave functions of u and in ideal mixing. The wave functions of the rj and rf 



It is also referred to as flavor-diagram or topological-diagram approach iiLthe literature. 
2 Historically, the quark-graph amplitudes T, C, E, A, P named in [10j were originally denoted by 
A, B, C, D, E, respectively, in 0, 0]. For the analysis of charmless B decays, one adds the variants 
of the penguin diagram such as the electroweak penguin and the penguin annihilation. 
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are given by 

\rj J \sm<p cos <p ) \ t] s ) 
where <j) = + arctan\/2, and 9 is the r]—rf mixing angle in the octet-singlet basis 

/ r/ \ _ / cos -sin0\/?78\ 
y r/' J \ sin cos 9 J \Vo J 

The most recent experimental determination of the rj — rj' mixing angle is eft = (40.4 ± 0.6)° from 
KLOE 



211 ] . 3 which is indeed close to the original theoretical and phenomenological estimates 



of 42.2° and (39.3 ± 1-0)°, respectively, made by Feldmann, Kroll and Stech [23j]. Notice that 
charm decays involving a final-state r/ s are governed by only one type of topological amplitudes, in 
addition to the associated hairpin diagram. For example, 

A{D° -> K°r) a ) = E, A{D° -> i^V) = E v , 

The D — > Mr] and Mr/ amplitudes then have the expressions 

A(D — )■ Mr/) = „4(L> -> Mr/ 9 ) cos <j> - A{D -> Mr/ S ) sin 0, 
«4(D->-Mr/) = «4(D -)• Mr/ g ) sin + -4(£> -)• Mr/ S ) cos c/>. (5) 

B. D ^ PP 

Based on flavor SU(3) symmetry, the quark-graph amplitudes for the Cabibbo-favored 
(CF), singly-Cabibbo-suppressed (SCS) and doubly-Cabibbo-suppressed (DCS) decays of charmed 
mesons into two pseudoscalar mesons are listed in Tables HI to IIIH respectively. Experimental data 
on the branching fractions are taken from the most recent CLEO measurements [4] . We shall follow 
the conventional practice to denote the primed amplitudes for SCS modes and double-primed am- 
plitudes for DCS decays. In the SU(3) limit, primed and unprimed amplitudes should be the same. 
Note that the selection rule for a vanishing Df — > tt + tt° follows from the isospin transformation 
properties of the weak Hamiltonian and isospin invariance of strong interactions and hence it is 



unaffected by SU(3) breaking or final-state interactions [241 ]. 

For D — > PP decays, the amplitudes T, C, E, A have dimensions of energy as they are related 
to the partial decay rate via 

Pc 



For the mixing angle (f> = arcsin(l/-\/3) = 35.26°, the rj and rj' wave functions have simple expressions 
22J: r] = ^(a/2?7 9 - Vs) = + dd- ss) and rj' = -^{Vg + V^Vs) = + dd + 2ss). 
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with p c being the center-of-mass momentum of either meson in the final state. The reduced quark- 
graph amplitudes T, C, E, A are extracted from the Cabibbo-allowed D — > PP decays to be (in 
units of 1(T 6 GeV) 

T = 3.14 ± 0.06, C = (2.61 ± 0.08) g^ 152 * 1 ) , 

E = (1.53±g;{g) e *( 122±2 )° , A = (0.39±g;S|) e^ 31+ -^° (7) 

for <f> = 40.4°, and 

T = 3.08 ± 0.06 , C = (2.46±°$) e - l(152±1) ° , 

E = (1.66 ± 0.06) e^ 120± V° , A = (0.34lg;g) e^ 70 ^" (8) 



for </> = 35.26°. The fitted x 2 value is 0.29 per degree of freedom with quality 59.2% and 1.73 
per degree of freedom with quality 18.8% for the fits flTJ and flSJ), respectively. Throughout this 
paper, we take A = 0.2253 ± 0.0007, the weighted average of CKMfitter 25J] and UTfit [26J groups, 
and neglect its small error. The relevant CKM matrix elements are V u d = V cs = 1 — A 2 /2 and 
V us = —V c d = A. With the inclusion of CKM factor V* s V u d in the amplitudes, our results in flSJ) 
agree with those in [18]. Note that the previous fit in 2008 had [itJ 

T = 2.78 ± 0.13, C = (2.04 ± 0.17) e *(- 151 ± 2 )° , 

E = (1.68 ± 0.12) e J ( 117±4 )° , A = (0.54 ± 0.37) e l ^ + -^° (9) 

with x 2 = 0.65 per degree of freedom. 

A few words about error extraction in our analysis are in order here. To obtain the 1-a range 
of some parameter around its best-fitted value, we scan the parameter while allowing all the other 
parameters to vary at the same time to minimize the x 2 function. This renders a more conservative 
range than the method that fixes all the other parameters to their best-fitted values when scanning 
the parameter of interest, as used, for example, in Ref. [ijj]. 

The magnitudes of the updated T and C are larger than the 2008 ones due to the newly measured 
branching fractions of K rj, it + r] and tt + t]' modes being larger than the old ones. Notice that the 
VF-annihilation amplitude A is substantially smaller than the VF-exchange one E, and it is no 
longer antiparallel to E, contrary to the previous claim ll|, Indeed, A is almost perpendicular 
to E in fit ([7]) and its angle relative to E is close to 50° in fit ([S]). 

A global fit to the data of singly Cabibbo-suppressed decay modes shows that the value of x 2 
is very large, of order 87 per degree of freedom. One solution we find is 

I* = 1.14, C = 2.36e i222 °, 

E' = 1.85e- i52 °, A' = 2.51e mo °. (10) 

This deviates substantially from the unprimed solution given in Eq. ([7]) . The main processes 
contributing to large x 2 are D° -> K + K~ (Ax 2 = 462), D° -)• vr+vr" (A^ 2 = 248), D° -)• 7r°7r° 
(Ax 2 = 119), and D+ (A X 2 = 96). This is an indication of large SU(3) breaking effects in 
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TABLE I: Branching fractions and invariant amplitudes for Cabibbo-favored decays of charmed 
mesons to two pseudoscalar mesons. Data are taken from Predictions based on our best-fitted 
results in ([7]) are given in the last column. 



Meson 


Mode 


Representation 


E 

(%) 


Hat 

(%) 


D° 


K-TT+ 

K°n° 
K° v 
K r\ 


V* s V ud {T + E) 

^V c * s V ud (C - E) 

VZV ud (C + E) cos - E sin </> ] 

V* s V ud [j- (C + E) sin + E cos <f> } 


3.91 ± 0.08 
2.38 ± 0.09 
0.96 ± 0.06 
1.90 ±0.11 


3.91 ±0.17 
2.36 ±0.08 
0.98 ±0.05 
1.91 ±0.09 


D + 


K°TT+ 


V* s V ud (T + C) 


3.07 ±0.10 


3.08 ± 0.36 


Df 


K°K+ 

vr+vr 

■K + r] 


V* s V ud {C + A) 


V* s V ud (V2A cos </> - T sin 0) 
V* s V ud (V2A smcj) + T cos 0) 


2.98 ±0.17 

< 0.037 
1.84 ±0.15 
3.95 ± 0.34 


2.97 ±0.32 


1.82 ±0.32 
3.82 ±0.36 



these decay modes. However, if we perform a fit restricted to the measured branching fractions of 
D° — > tt°7], irrf , f]7], rjrf , the extracted amplitudes are (in units of 10~ 6 GeV and for <fi = 40.4°) 

C' = 2.88±g;if, E ' = (l-37±g;i|) e -< 79±5 )° , (11) 

with x 2 = 0.29, where we have assumed C to be real and the phase on E is relative to C. Note 
that there is a two-fold ambiguity in the relative strong phase, and we have picked the negative 
one that is consistent with Eq. (J7J) . These values agree with the unprimed values in Eq. ((7J) within 
l-o". If the data of D — > ttV is also included in the above fit, it will dominate the x 2 value which 
now reads 19.6. A global fit to both CF and SCS modes yields x 2 = 52 per degree of freedom. 

Assuming that the primed amplitudes are identical to unprimed ones, i.e., T' = T, C = C, 
etc., the predicted branching fractions of SCS D —> PP decays are shown in the last column of 
Table |TTJ It is clear that the predicted rates for -k + it~ and 7r°7r° are too large, while those for 
K + K~, 7T + 7r°, iT + rj, tt + t]', K + rj and K + rf are too small compared to experiments. The decay 
D° — > K°K° is prohibited by SU(3) symmetry, but the measured rate is comparable to that of 
D° —7- 7r°7r°. In the next section we shall study flavor SU(3) symmetry violation in SCS decays 
using the factorization approach and investigate its possible connection to final-state interactions. 
Setting the double-primed topological amplitudes as unprimed ones, the branching fractions of 
DCS processes are displayed in Table IIII1 

A few remarks on our assumption of ignoring the SE and SA amplitudes are in place here. The 
former (latter) diagram appears only in the singly Cabibbo-suppressed decay modes D° — > 7T ?/') 
and D° — > rp]^ (D + — > Tr + r]^ and Df — > K + r]^). Ref. [l8|] includes them in their analysis. 
Two solutions have been found for the magnitude of SE. One is relatively large when compared 
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TABLE II: Branching fractions and invariant amplitudes for singly-Cabibbo-suppressed decays of 
charmed mesons to two pseudoscalar mesons. Data are taken from Q]. Predictions based on our 
best-fitted results in (|7|) with exact flavor SU(3) symmetry are given in the last column. 



Meson Mode Representation 


B 


^theory 




(xl(T 3 ) 


(xl(T 3 



D° 



V c * d V ud (T' + E>) 

j= 2 v; d v ud (c - e>] 

-V^V ud E'co S cj>- 



^V* s V us C'sm ( 



-V* d V ud E' sin <j> + ^V* S V US C cos , 



7T 1 7T 
7T°7T 

TT°r] 

r)V -^V c * d V ud (C + E>) cos 2 </, + V* S V US (2E' sin 

VV' -lV^V ud (C + E')sm2(f) + V c * s V us (E'sm2<j)- -^C" cos 20) 1.05 ±0.26 1.19 ±0.07 

K+K- V* S V US (T' + E') 4.07 ±0.10 1.92 ± 0.08 

K°K° V* d V ud E' s + V* s V us E' d a 0.64 ±0.08 



1.45 ±0.05 2.24 ±0.10 

0.81 ±0.05 1.35 ±0.05 

0.68 ±0.07 0.75 ±0.02 

0.91 ±0.13 0.74 ±0.02 

^C" sin 20) 1.67 ±0.18 1.44 ±0.08 



* + *° ^V c Wud{T> + C>) 

-^V* d V ud (T' + C + 2A>) cos 4> - V*V US C sin , 
j=V c * d V ud (T' + C + 2A>) sin + V*V US C> cos . 
V c * d V ud A' + V c * s V us T' 



ir + rj' 



1.18 ±0.07 0.88 ±0.10 
3.54 ±0.21 1.48 ±0.26 
4.68 ±0.30 3.70 ±0.37 
6.12 ±0.22 h 5.46 ±0.53 



D+ tt + K° V* d V ud T' + V* S V US A' 2.52 ± 0.27 c 2.73 ± 0.26 

tt°K+ -^{V* d V ud C -V* S V US A') 0.62 ±0.23 0.86 ± 0.09 

if +77 -t=(V* d V ud C' + V^V us A')cos<p-V* s V us (T' + C + A')sin0 1.76 ±0.36 0.78 ± 0.09 

K+r]' ^(V* d V ud C' + V* s V us A')sm(j) + V* s V us (T' + C' + A')cos0 1.80 ±0.52 1.07 ±0.17 

°The subscript q in E' q refers to the quark-antiquark pair popping out of the-wacuum in the final state. 
b A new Belle measurement yields B{D+ K+K°) = (5.50 ± 0.16) x lO" 3 ^. 



C A new Belle measurement yields B(D+ tt+K°) = (2.40 ± 0.18) x 10" 3 [27|. 



to the other amplitudes (T, C, E, and A) and thus argued unlikely. The other is close to zero, 
and explains our good fit of C and E' given in Eq. (jlip even without including SE. As to the 
magnitude of SA, Ref. 3] finds only one large solution that is again unlikely in comparison with 
the other amplitudes. Moreover, our predictions for the branching fractions of the D + — > tt + i]^ 
and Df — > Krj-'^ decays in Table Hi] which are based on the fit to Cabibbo-favored modes without 
including the SA amplitude, are in better agreement with the data than the predictions given in 
Table VI of Ref. [18| where the large SA solution is used. Therefore, the current data do not 
suggest significant SE and SA contributions. 

The quantities R(D°), R{D+) and R(D+) defined by 

T(D° K s ir°) - T(D° K L ir°) 



R(D ) 



T{D° -> K s tt°) + T{D° -> K L 7T°) 
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TABLE III: Branching fractions and invariant amplitudes for doubly-Cabibbo-suppressed decays 
of charmed mesons to two pseudoscalar mesons. Data are taken from Predictions based on our 
best-fitted results in (|7|) with exact flavor SU(3) symmetry are given in the last column. 



Meson 


Mode 


Representation 




l xiU 1 


^theory 
{X 1U ) 




II /I 

K°7T° 


V* V (T" -A- F,"\ 
-±V*V US (C" - E") 




x.^io zn u.u t 


1 i o _|_ n nc: 

0.67 ±0.02 




K°r] 


V^V us [^(C" + E")co S ^ 


b - E"sinqb] 




0.28 ±0.02 




K°rj' 


V c * d V us [j= 2 (C» + E")sm q 


> + £"cos<£] 




0.55 ±0.03 


D+ 


K°TT + 


V c * d V us (C" + A") 






1.98 ±0.22 




K+TT° 


j= 2 V c * d Vus(T" - A") 




1.72 ±0.19 


1.59 ±0.15 




K+T] 


V: d V us (^(T" + A")cosc 


b - A" sin 4>) 




0.98 ±0.04 




K+7]' 


V^V us (^(T" + A») S m q 


) + A" cos 4>) 




0.91 ±0.17 


Df 


K°K+ 


V c * d V us (T" + C") 






0.38 ± 0.04 



+ T(D+ -» K 5 vr+) - Tp+ -> if L 7r+) 

1 J r(z?+ -> i^ 5 7r+) + r(L>+ ' 
Rm+ ^ _ r(g+ -> - rQD+ -» k l k+) 

T(Df^K s K+) + T(Df^K L K+) ^ 
measure the asymmetries between K$ and i?x production in these decays due to the interference 
between CF and DCS amplitudes. It is straightforward to show that 4 

R(D°) = 2 tan 2 9 C = 0.107, 

/ c" + A" \ 

R(D + ) = 2 tan 2 9 C Re f - - - j = -0.019 ± 0.016 , 

/ n" + T" \ 

R(D+) = 2tan 2 6> c Re( - ^ - j = -0.008 ± 0.007 . (13) 
Note that R(D + ) is predicted to be 0.035 ~ 0.044 in [29] and -0.005 ± 0.013 in 0, while R(D+) 



is -0.0022 ± 0.0087 in [lg]. Our results for R(D + ) and R(D+) differ in central values from those 
in [ijj since the latter were obtained using the topological amplitudes extracted for <p = 43°. The 



The formula R(D°) = 2 tan 2 8c was originally derived under the phase convention such that Ks 



-^(K — K°) and Kl = + K°) [28|- Of course, physics should be independent of the choice of 

phase convention. If we choose the convention CP\K°(p)) = \K°(—p)), we will have Ks — + 
K°) and Kl = -j=(K° — K°) in the absence of CP violation. Since the axial vector current = 
S7 M (1 — j 5 )d transforms as CP(s7 M (l — j 5 )d)(CPy = —drf^l — j 5 )s under the CP transformation, it 
is clear that the decay constants of K° and K° are opposite in sign. Consequently, A(D° — > K°ir°) — 
-(V£tV us )/(V* s V ud )A(D° -> K°ir°) = tan 2 9 C A(D° -> K tt°). It follows that D° -> K°ir° and D° -> 
K°tt° contribute constructively to D° — > Kgir and destructively to D° — > K^ir . 
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CLEO measurements are 30] 



R(D°) = 0.108 ±0.025 ±0.024, R(D + ) = 0.022 ± 0.016 ± 0.018. (14) 

Evidently, the calculated R(D ) agrees well with experiment, while there is no evidence for a 
significant asymmetry in D + — > Kg^^ decays. 



C. D^VP 

The topological amplitude expressions for CF, SCS and DCS D — > VP decays are listed in Tables 
IIVI to IVll respectively. For reduced amplitudes T and C in D — > VP decays, the subscript P (V) 
implies a pseudoscalar (vector) meson which contains the spectator quark of the charmed meson. 
For E and A amplitudes with the final state qiq~2, the subscript P (V) denotes a pseudoscalar 
(vector) meson which contains the antiquark q<i . 

The invariant amplitude is related to the partial width via 

r ^ yp i = 4»f < 15 > 

u pol. 

by summing over the the polarization states of the vector meson, or through the relation 

by taking the polarization vector out of the amplitude, where A = (my /mo)A (e ■ po)- The first 



approach has been used in 36|] for the extraction of topological amplitudes in B — > VP decays, 
while the second approach was employed in [ll|, 3] for D — > VP . To illustrate the relations 
between Eqs. (|15p and (|16p . we consider the channel — > tt + 4> as an example. Its amplitude is 
given by A = V* s V u dTy with 

Ty = 2UmtA$'+{ml)(e ■ p D ). (17) 

evaluated in the factorization approach. After summing over the polarizations of the vector meson, 
the above expression can be simplified by replacing my e-pu with moPc- Hence, we obtain Eq. (|16() 
with 

f v = 2Um D A^(ml) = ^^~- (18) 

my £ ■ po 

From the decays Df — > tt + 4>, D° — > K°(f> and the three CF D — > K*ir channels, two solutions for 
the magnitudes of and relative phases between Ty, Cp and Ep are found and shown in Table rvTIl 
for two different ways of extracting topological amplitudes using (|15|) and (|16|) . It is noted that 
these results have no dependence on the 77-7/ mixing angle. Solutions with a smaller Cp, i.e. A' 
and S', are ruled out by the measurements of SCS decays D° — > ir°(j) and D + — > ir + (fi since they 
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TABLE IV: Branching fractions and invariant amplitudes for Cabibbo-favored decays of charmed 
mesons to one pseudoscalar and one vector meson. Data are taken from [3l|. The last two columns 
show the fitted branching fractions obtained from solutions (A, Al) and (S, SI), respectively, for 
<j> = 40.4°. Due to the lack of information on Ap and Ay, no prediction is attempted for Df decays 
except Df — > ir + (f). 



Meson Mode 


Representation 




'-'cxp 


Bfit (A,A1) 


Bat (S,S1) 










(%) 


(%) 


(%) 


D° 


K*-TT + 


V* s V ud {T v + E P ) 




5.91 ±0.39 


5.91 ± 0.70 


5.91 ±0.66 




R-p + 


V*V ud (T P + E v ) 




10.8 ±0.7 


10.8 ±2.2 


10.7 ±2.3 




K V 


^V* s V ud {Cp - Ep) 




2.82 ± 0.35 


2.82 ±0.34 


2.82 ±0.28 




K p° 


■^V* s V ud (C v - E v ) 




1.54 ±0.12 


1.54 ± 1.15 


1.55 ±0.34 




— *0 

K rj 


V*V ud (^(C P + E P ) cosq 


*) — Ev sin <j 


i) 0.96 ±0.30 


0.96 ± 0.32 


1.12 ±0.26 




k V 


V*V ud {±(C P + Ep)sm4 


) — Ey COS (, 


J) <0.11 


0.012 ± 0.003 


0.020 ± 0.003 




-rrO 

K u 


±V c * s V ud (Cv + E V ) 




2.26 ±0.40 


2.26 ± 1.38 


2.34 ±0.41 




i\ Cf> 


V c * s V ud E P 




0.868 ±0.060 0.868 ±0.139 


0.868 ±0.110 


D+ 


17*° + 


V* s V ud (T v + C P ) 




1.83 ±0.14 


1.83 ± 0.49 


1.83 ± 0.46 




K p* 


V* s V ud (T P + C v ) 




9.2 ±2.0 


9.2 ±6.7 


9.7 ±5.2 


Df 


K*°K+ 
K°K*+ 

p + 7T° 


V c * s V ud (C P + A v ) 
V* s V ud (C v + A P ) 
-t=V c * s V ud (A P - A v ) 




3.91 ±0.23 a 
5.3 ± 1.2 








p + 7] 


V* s V ud (^(Ap + A v )cos< 


b — Tp sin 4 


>) 8.9 ±0.8 6 








p+rj' 


V*V ud (^(A P + A v )smq 


) + Tp cos 4 


>) 12.2 ±2.0 








vr+p° 


j=V c * s V ud (A v - A P ) 














±V* s V ud (A v + A P ) 




0.21 ±0.09 c 








TT + (f) 


V* s V ud T v 




4.38 ±0.35 


4.38 ± 0.35 


4.38 ±0.35 



Trom the fit fraction T(D^ 
PDG value, [m |- 



6 Prom laa. 

c From 



K*°K+)/T(D+ -> K+K- 



(47 .4 ± 



1.5)% H 



combined with the 



will lead to predictions too small by a factor of 3 when confronted with experiment, 5 while the 
data can be nicely explained with the larger Cp (see Table W\ , as noticed in [l3.[l3|. Although the 



5 There are two different values of B(D° -> 7r°(/>) quoted by the Particle Data Group |3lJ] : (1.24±0.12) x 10 



-3 



obtained by BaBar [37] and CLEO 38] from the Dalitz-plot analysis with interference, and (0.76±0.05) x 
obtained by Belle by measuring the background for the radiative decay D° — > <jyy [39j. In this work, 
we shall use the former. 
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TABLE V: Same as Table IIVI except for singly-Cabibbo-suppressed decays of charmed mesons. 
Due to the lack of information on Ap and Ay, no prediction is attempted for Df decays except 
D~ 



7T 



Meson Mode 


Representation 


B 

^exp 


^theory (A,Al) ^theory (S,S1) 








(xlO" 3 ) 


(xl0~ 3 ) 


(xl0~ 3 ) 


D° 


7T + p 


V c * d V ud (T v + E> p ) 


4.97 ±0.23 


3.92 ± 0.46 


5.18 ±0.58 




1 

7T /r 


V c * d V ud (T' P + E' v ) 


9.8 ±0.4 


8.34 ± 1.69 


8.27 ± 1.79 




TT°p° 


2 V cdYud{C'p + C'y ~ E'p - E'y) 


3.73 ± 0.22 


2.96 ±0.98 


3.34 ±0.33 




K K 


V*V US (T V + E'p) 


1.53 ±0.15 


1 on i n o A 
1.99 ± 0.24 


1 nn _i_ n on 
1.99 ± 0.22 




If — 
1\ I\ 


V CS V US {T P + Ey) 


4.41 ± 0.21 


A 9^ -I- n Sfi 
4. ZD It U.oD 


o.io it u.oy 




K°K*° 


V CS V US E P + V cd V ud E v 


< 1.8 


0.29 ±0.22 


0.05 ± 0.06 




i\ i\ 


V c * s V us E'y + V c * d V ud E' P 


< 0.9 


n 9Q -i- n 99 
u.zy m u.zz 


U.UO nz U.UD 






~2 V cd V ud{^y - Op + iip + tjy) 


^ n or 

<. U.ZD 


0.10 ±0.18 


1.01 ±0.18 




TT°(j) 


-J^VcsVusC'p 


1.24 ±0.12 


1.22 ± 0.08 


1.11 ± 0.05 




rjuj 


"KsK^C^sin^ 


o 91 _|_ n oq a 


3.08 ± 1.42 


3.94 ±0.61 




■q'uj 


^K4(C(, + C'p + E' v + E'p) sin 

+ KsK S ^C( / COS,/» 




0.07 ±0.02 


0.15 ±0.01 




T]4> 


V*V US (^C' P cos - + E'p) sin 


) 0.14 ±0.05 


0.31 ±0.10 


0.41 ±0.08 




7]p° 


KdK4(Cv -C'p-E'y- E'p) cos 
-K^Ks^C^sin^ 




1.11 ±0.86 


1.17 ±0.34 




ri'p 


V* d V ud \(C' v -C'p-E'y- E'p) sin 




0.14 ±0.02 


0.26 ±0.02 



D+ 


TT+p° 
7T°p+ 


j= 2 V; d V ud (T v + C'p-A' P + A'y) 

j= 2 V c * d V ud (T'p + C'y+A'p-A'y) 


0.82 ±0.15 




it + lj 


^= 2 V c * d V ud (Ty+C'p + A'p + A'y) 


2.1 ±0.9 6 




ir + (f) 


V* s VusCp 


6.2 ±0.7 6.21 ±0.43 5.68 ± 0.28 




VP + 


VcWud^T'p + C'y + A'y + A'p) cos <j> 

-V*V US C V sincp 


< 7 




rj'p + 


y^ud^T'p + C' V + A'y + A'p) sin0 


< 5 




K+K*° 


+V*V US C V cos (P 






V c * d V ud A'y+V c * s V us Ty 


4.4 ±0.5 




K°K*+ 


V c * d V ud A'p + V c %V us T'p 


31.8 ± 13.8 


Df 


ir+K*° 


V* d VudT v + V* S V US A' V 


2.25 ±0.39 



7T°K* + 

K+p° 
K°p+ 

T]K* + 
7]'K* + 

K + (p 



'cs Vus-A-y 

V c * d V ud T'p + V c * s V us A' P 



^^udC v -v; s v us A 
j= 2 (V: d VudC'p-V* s V us A 



'p) 



2.7 ±0.5 



v cd v ud- L p^~ v cs y us^p 

j= 2 (V; d V ud C v + V* S V US A' V ) cos 

-V* S V US (T P + Cy + A'p) sine 
-^(V: d V ud C v + V C * S V US A'y) sin 

+V^V US (T' P + C v + A' P ) cos, 
^(V; d V ud C' P + V*V us A P ) 

V* s V us (Ty+C'p+A'y) 

12 



< 2.4 h 

< 0.57 



"Data from 



35 1 
4/11 



TABLE VI: Same as Table ITVl except for doubly-Cabibbo-suppressed decays of charmed mesons. 



Meson Mode Representation S exp ^theory (A,A1) ^theory (S,S1) 

(xlQ- 4 ) (xlQ- 4 ) (xlQ- 4 ) 

D° K* + tt- V* d V ua (T£ + E' v ) 3.0^| 3.59 ±0.72 2.69 ± 0.58 

K*°tt° -±=V* d V us (Cp - E v ) 0.54 ±0.18 0.74 ±0.17 

<pK° V* d V us E' v 0.06 ±0.05 0.15 ±0.06 

p~ K+ V* d V us (T v + E P ) 1.45 ±0.17 1.91 ± 0.21 

p°K° ^V* d V us (C v - E' P ) 0.91 ±0.51 0.63 ±0.19 

coK° -^V* d V us (C'{r + E' P ) 0.58 ±0.40 0.85 ± 0.21 

K*°r] V* d V us (^{C' P + E' V ) cos <j) - E' P sin (f>) 0.33 ± 0.08 0.28 ± 0.05 

K*° Tj' V* d V us ( 75 (Cp + E v) sin (P + e p cos </>) 0.0040 ± 0.0006 0.0061 ± 0.000 4 

D+ K*°tt + V* d V us (C P + A v ) 4.35 ±0.90 

K* + 7T° j=V* d V us (Tp'-A>{ / ) 
<PK+ V* d V us A' v 
p+K° V c * d V us (C v + A' P ) 

p o K+ -^V c * d V us (T» - Ap) 2.4 ±0.6 

j= 2 V: d V us (T»+Ap) 
K* + V V* d V us (-^(Tp + Ay ) cos 4> - A" P sin 0) 
K* + rl V* d V us (MTp- + A v ) sin cf> + A" P cos </>) 



D+ K* + K° 


V; d V us (T' P + Cy) 


1.17 ±0.86 


1.03 ±0.55 


K*° K+ 




0.20 ±0.05 


0.22 ±0.06 



fact that the magnitude of color-suppressed Cp is larger than the color-allowed Ty seems to be in 
contradiction with the naive expectation, we shall see in the next section that the corresponding 
effective parameter a\ is still larger than a 2 , as anticipated from the short-distance approach. 

Using the solutions of Ty, Cp and Ep as inputs, the other amplitudes Tp, Cy and Ey can 
be obtained by considering the decay modes D° — > K°co, D° — > K*°r], D° — > K~p + , K°p° and 
D + —7- K°p + . Here we have assumed that Tp and Ty are relatively real. We obtain 6 solutions 
Al— A6 for topological amplitudes extracted from (|16p and 5 solutions SI— S5 when amplitudes 
are extracted from (|15p . as given in Table IVTTT1 

Several remarks are in order: (i) All the solutions are exact with \ 2 = 0, except for SI where 
the fit x 2 value is 0.42. This is reflected in some of the predicted branching fractions in Table HVl 
(ii) Solutions displayed in Table IVTHl are obtained using the r) — rf mixing angle <j) = 40.4°. If the 
angle (j) = 35.26° is employed, we obtain 

(Al) T P = 7.85t° 2 ii C v = (3.63±S;g) e^ 172+ ^° , Ey = (2.51^;^) e -( 11 o^)° , 

(SI) T P = 3.12j££, Cy = (1.251°$) e ^ 77 ^?)°, Ey = (l^l^e^ 107 ^) . (19) 

where the amplitudes in the Al solution are quoted in units of 10 -6 and those in the SI solution 
in units of 10 _6 (e • pd)- By comparing them with solutions Al and SI in Table lyTTlj we see that a 
small decrease of the 77—7/ mixing angle will reduce Tp and Cy slightly and enhance the magnitude 
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TABLE VII: Solutions for Ty, Cp and Ep. The top two solutions in units of 10 are obtained 



using Eq. (fT6|) . while the bottom two solutions in units of 10 e (e-pD) are extracted using Eq. (fT5]) . 





T v 






£7p 


A 


4 16+ 016 
4 - iD -0.17 




-i(162±3)° 


(3.09±0.11)e^ 93±5 )° 


A' 


4 1fi +ai6 

4 - 1D -0.17 


(3.00i°;l)e- 


-i(158±2)° 


(S.OOiO.ll)^^ 5 ) 


S 


o i rr+0.08 
z - id -0.09 


(2.68±tiit)e- 


i(164±3)° 


(1.69 ± 0.06)e- i ( 103±4 )° 


S' 


9 i rr+0.08 


(l-53i°i 5 7 )e- 


-i(162±j?)° 


(1.69±0.06)e i ( 103±4 )° 



TABLE VIII: Solutions for Tp, Cy and Ey with inputs from solution A and S in Table IVIIL 
Solutions Al— A6 in units of 10 -6 are obtained using Eq. (|16p . while solutions SI— S5 in units of 
10 _6 (e • pd) are extracted using Eq. (|15p . Here we take cp = 40.4°. 





T P 


Cy 


-By 




Al 


o -I -,+0.32 




(l-51^i)e- 


<(124tSS)° 


A2 


6.i6^;l 




(3.9Bj*g)e 


-i(89±6)° 


A3 


-6.5918;™ 


(4.09i°-l) e - i ( 39 -u)° 


(l-66±o^)e- 


i(113±16)° 


A4 


-6.39^:81 


(2.19t 1 - ?)e- i ( 12 --)° 


(3.83i?;^)e- 


-i(89±10)° 


A5 


o no+0.34 
°-°°-0.37 


(1.74^-^) e -( m -")° 


(4.06i°;l)e- 


i(185±^)° 


A6 


_o 4fi +0.39 

-j- 4D -o.45 


(2.03t^ 9 )e i ("-»)° 


(3.93l°; 3 «)e- 


i(185±10)° 



SI 

S2 
S3 
S4 
S5 



o 14+O.3I 

_ 9 9n+0.96 
^ u -0.35 
9 i 1 +0.64 
z - il -0.37 

_i 41 + - 17 
i - 41 -0.30 



-1.38 



+0.14 
-0.15 



(1.33±g;|?)e i ( 177 iii) 



(0.48t£g)e 
(0-82i°; 34 )e 
(0.75±g : |f)e 



-i(53llS)° 



+0.43^-i(39±gg)° 
i(102t^)° 



(1.31^- 4 ?)e- J ( 106 -i6)° 

(l.isig-De-^ 131 -") 

(1.79^i2) e -(i04t? 2 )° 
(1.66^-^ )e -(i77^)° 

(1.69±g-}t)e- < C 177 ±"^ 



of By. (hi) Using cp = 35.26° and removing the CKM angles V* d V u d from the amplitudes, we are 



able to reproduce solutions Al— A6 in Table III of 



3. 



(iv) The relation Ey = —Ep used in 



[l^ is not borne out in this analysis. Instead, the angle between Ey and Ep turns out to be 
small, (v) To see the relative phase St v t p between Tp and Ty, we fit to the data of the CF D°, 
D + decays and the channel — > 7r + cp and find that St v t p is of order (16 ~ 18)° when using 
Eq. (|15p and (—1 ~ 1)° when using Eq. (|16p . Therefore, the assumption of relatively real Tp and 
Ty amplitudes is justified, (vi) A careful comparison between the two sets of solutions (Al— A6 and 
S1-S5) in Table ELLD shows that there is no solution in the latter set that corresponds to solution 
A2. Ref. 16] claims that solution A2 is an alternative solution to Al as it produces roughly the 
same x 2 as Al when they are applied to the SCS decays. However, we find that the D° — > cprj 
mode alone contributes Ax 2 = 285 in solution A2, sufficient to disregard this solution. 



6 It seems to us that there is a sign typo in solution A6 of [16]: The amplitude Tp there should be negative. 
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Now we come to the remaining parameters Ap and Ay appearing in the CF decays. The 
non-observation of the decays Df — > and /? + 7T° gives a strong constraint on the TV-annihilation 



amplitudes. The mode Df — > p°ir + was quoted by the Particle Data Group as "not seen" [31]. 
While the ratio T(D+ -> p°ir + )/T{D+ -> vr+yr+vr-) was measured by E791 to be (5.8 ±2.3 ±3.7)% 
I40I ]. a most recent Dalitz plot analysis of Df — > ir + ir + ir~ by BaBar yielded the fit fraction 



(1.8 ± 0.5 ± 1.0)% [4l|]. This mode was also not seen by the FOCUS Collaboration [4g|. Given 



the branching fraction B(Df — > tt + tt + tt^) = (1.11 ± 0.08)% [3l|, the BaBar result leads to the 
upper limit B(Df — > p°ir + ) < 5 x 10~ 4 . If we ignore the decays — >• p°ir + and p + ir° for 
the moment, the annihilation amplitudes Ap and Ay can be determined from the four channels 
K*°K + , K°K* + , p + r) and 7r + w in conjunction with the information of Cp from solution A or S 
and Tp, Cy from solution Al or SI. It turns out that although the above four Df modes are 
nicely fitted, the predicted B(D+ — > of order 4 x 10 3 is too large. 

From the measured rates of Df — > 7r + p° and Df — > tt + oo we have 

\A V - A P \ < 0.20 x 10~ 6 (e -p D ), \A V + A P \ = (0.41 ± 0.09) x 10~ 6 (e • p D ) (20) 

obtained using Eq. (TT5]) and 

\A V -A P \ < 0.50 x 10~ 6 , \A V + A P \ = (1.04 ±0.22) x 10~ 6 (21) 

using Eq. (fl~6j) . Hence, the magnitudes of and Ay are much smaller than that of Cp and Cy. 
In principle, the annihilation amplitudes can be determined when the above equation is combined 
with the amplitudes of Df -)• K*°K + and K°K* + . However since \C P \ > \Cy\ > \A P \, \A V \, it 
is not possible to have a nice fit to the data of K*°K + , K°K* + , tt + lo and 7r + p° simultaneously. 
To see this, we set Ay ~ Ap as a consequence of the non-observation of Df — > Under this 

relation, it is naively expected that Df — > K®K* + has a rate larger than — > K*°K + since 
\C P \ > \C V \. Experimentally, it is the other way around: B(D+ K*°K + ) = (3.97 ± 0.21)% 
and B(Df — > K°K* + ) = (5.3 ± 1.2)%. Since the former has been measured several times in the 



past decade (see Particle Data Group [3l|] and the latest one [32|]) while the latter was measured 
two decades ago [4^], it is likely that the quoted experimental result for Df — > K°K* + was 
overestimated. In short, a sensible determination of Ap and Ay cannot be made at present and 
we have to await more accurate measurement of — > K°K* + . 

The cited experimental result for Df — > p + rj' is also problematic. Since \Tp\ ^> \ Ay ± Ap|, we 
can neglect the annihilation contributions for the moment. It follows that if we use the solution with 
the largest T P , the predicted branching ratios will be B{D+ -)• p+r?) = (7.4 ± 0.7)% [(7.1 ± 1.8)%] 



7 In Ref. [16|, the amplitudes A v and A P were determined from D+ (K*°K + , K°K* + . tt+uj). It is 
not clear to us how the four unknown parameters (two magnitudes and two pha ses) can be extracted 
out of three data points. Anyway, solutions Al and A2 given in Table IV of 16( for Ap and Ay lead 
to B(D+ p°tt + ) ~ 1.2% and 2.2% in their central values, respectively, which are obviously too large 
compared to experiment. 
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and B(D+ p+rf) = (2.7±0.2)% [(2.6±0.7)%)] for solution Al [SI]. While the predicted rate for 
the former is close to the recent CLEO measurement B(D^ — > p + rf) = (8.9 ±0.8)%, the latter is far 
below the quoted result B{Df — > p + r]') = (12.2 ± 2.0)%. It is very unlikely that the flavor-singlet 
contribution unique to the r/o production can enhance its branching fraction from 3% to 12%. This 
issue should be clarified by new measurement of Df — > p + r/' . 

Among the solutions Al— A6 and SI— S5, Al and SI are our preferred solutions for the following 
two reasons: (i) A global fit to the SCS decays of the D° meson (data of D + and Df are not included 
in the fit due to the lack of information on Ap and Ay) indicates that Al and SI have the lowest 
values of x 2 > 89 and 165, respectively, and (ii) Tp has to be sufficiently large in order to possibly 
accommodate the data of — > p + rj. 



III. PHENOMENOLOGICAL IMPLICATIONS 



A. Parameters a\ and a 2 

In the diagrammatic approach, the topological amplitudes extracted using the x 2 fit method 
are not unique. Indeed, several other solutions are also allowed. It is useful to consider the 
amplitudes for color-allowed and color-suppressed diagrams in the factorization approach in order 
to discriminate between different solutions. For CF D — > Kit, K*tt, Kp decays, the factorizable 
amplitudes read 



^(Kn) U{ml - m 2 K )F DK (m 2 



C = ^a 2 (Kn)f K (ml-ml)F D *(m 2 K ), 
T V = ^a 1 (K*ir)2f 7r m K *A{? K *(ml)(e-p D ), 

C P = ^a 2 (K\)2f K *m K *F 1 D7r (m 2 K ,)(e- PD ), (22) 



^ 1 (Kp)2f p m p F 1 1K iwin:-pn) 



Gp 

Cy = -j=a 2 {Kp)2f K m p A^ p {m 2 K ){£ -p D ), 

where we have followed [44] for the definition of form factors. Factorization implies a universal, 
process-independent oi and a 2 , for example, a 2 (Kp) = a 2 (K*ir) = a 2 (Kir). 

Under the naive factorization hypothesis, one has 

1 1 

ai(/x) = ci(p) + — c 2 (/x), a 2 (p) = c 2 (p) + —c^p), (23) 

for decay amplitudes induced by current-current operators Oi t2 (p), where c\ )2 {p) are the corre- 
sponding Wilson coefficients and N c is the number of colors. However, it is well known that this 
naive factorization approach encounters two difficulties: (i) the coefficients aj given by Eq. (|23p are 
renormalization scale and 75-scheme dependent, and (ii) it fails to describe the color-suppressed 
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TABLE IX: Form factors at q 2 = and the shape parameter a. 







tpDK 


pD-K 


pDK 


A?" 


A™* 


F(Q) 
a 


0.666 
0.21 


0.739 
0.30 


0.666 
0.24 


0.739 
0.33 


0.74 
0.36 


0.78 
0.24 



decay modes due to the smallness of a 2 . In particular, the ratio of D° — > K°ir to D° — > K~ir + is 
predicted to be of order 10 -2 , while experimentally it is close to 1/2. Therefore, it is necessary to 
take into account nonfactorizable corrections: 



ax 



c 1 (/i) + c 2 (M) ( jjr+Xi(M) 



1 



u-2 ~ <-o(//) + 'Ml//) ( jy - +X2(M) 



(24) 



where nonfactorizable terms are characterized by the parameters Xi> which receive corrections 
including vertex corrections, hard spectator interactions involving the spectator quark of the heavy 
meson, and FSI effects from inelastic rescattering, resonance effects, etc. The nonfactorizable terms 
Xi(p) will compensate the scale and scheme dependence of Wilson coefficients to render a% physical. 
In the so-called large- N c approach, a rule of discarding subleading 1/N C terms is formulated Q]. 
This amounts to having a universal nonfactorizable term xi = X2 = — 1/N C in Eq. (|24p and hence 

d(m c ) = 1.274, a 2 » c 2 (m c ) = -0.529 (25) 



for A^g = 215 MeV and m c {m c ) = 1.3 GeV [45] . This corresponds to a relative strong phase of 
180°. Empirically, this set of a\ and a 2 gives a good description of the hadronic charm decays. 
Hence, a\ and a 2 in the large- N c approach can be considered as the theoretical benchmark values. 

From the topological amplitudes obtained in Sec. II, we are ready to determine the effective 
Wilson coefficients a\ and a 2 and their ratios. For the invariant amplitudes of D — > VP determined 
from Eq. (|15p . a\^ 2 are extracted using the factorizable amplitudes given in Eq. ()22|) . However, if 
the topological amplitudes of D — > VP are extracted from Eq. (|16fl . we should use the following 
factorizable amplitudes 



C P 
f P 
Cv 



ax{K ir)2f 7r m D A% 



DK* 



Op 

V2 
Gp 

72 

^a 1 (Kp)2f p m D F 1 DK (m 2 p ) 



a 2 {K^)2f K *m D F^(m A K *), 



[m 2 K ) 



(26) 



^=a 2 (Kp) 2f K m D AQ P (i 

to determine 01,2- As noted in passing, the amplitudes A and A are related by Eq. (fT8|) . 

In order to extract the parameters ai j2 we need to know the form factors and their q 2 dependence. 
There exist many model and lattice calculations for D to n, K transition form factors. In this work 
we shall use the following parametrization for form-factor q 2 dependence 46J] 

F(0) 



(1 - q 2 /m 2 )(l - aq 2 /m 2 ) 



(27) 
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TABLE X: The extracted parameters a\ and a<i- 







D - 


+ k\ 


D - 




A 


S 


Al 


SI 


l°ll 


1.22 ±0.02 


1.32 ±0.05 


1.43 ±0.06 


1.38 ±0.06 


1.29 ±0.17 


\ a 1\ 


0.82 ±0.02 


0.87 ±0.05 


0.95 ±0.05 


1.06 ±0.12 


0.82 ±0.27 


1 02/01 1 


0.67 ±0.02 


0.66 ± 0.05 


0.67 ±0.05 


0.77 ±0.09 


0.63 ±0.22 


Arg(a 2 /ai) 


-(152 ± 1)° 


-(162 ±3)° 


-(164 ±3)° 


(164 ± 23)° 


(177 ± 15)° 



with m* being a pole mass. Specifically, m* = mo*, mc,, mo*, for form factors F®^ , Aq K * , 
and Aq P ', respectively. The inputs for the parameters F(0) and a are summarized in Table 
IIXI Form factors for D to tt and K transitions are taken from the recent CLEO-c measurements 
of D meson semileptonic decays to tt and K mesons [47] , while D — )• p, K* transition form factors 
:'rom [46|] with some modification. For decay constants we use f n = 130.4 MeV, Jk = 155.5 MeV 



3l|, f p = 216 MeV and f K * = 220 MeV [4J]. 

The extracted parameters are shown in Table [Xj where we have used solutions A and S 
for the topological amplitudes Tp,Cy and Tp, Cy, respectively, and the preferred solution Al for 
Ty,Cp and SI for Ty, Cp. Due to final-state interactions, the magnitude of a 2 is substantially 
larger than the benchmark value 0.529 , while the magnitude of a\ is of order 1.3 ~ 1.4 for D — > K*tt 
and D — > Kp decays. 

As noticed in passing, solutions of Cp with magnitude smaller than Ty (solutions A' and S' 
in Table |VlI]) are ruled out since the predicted rates for SCS decays D° — > ir (f) and D + — > -K + (j) 
are too small compared to experiments. This can also be seen from the extracted ratio |«2/oi| 
which is expected to be similar for all two-body hadronic D decays. For solutions A' and S' with 
smaller Cp, the ratio (02/01! is 0.38 and deviates substantially from the value of 0.68 obtained 
from D — > Kir decays. Notice that although \Cp\ > \Ty\, the extracted parameters look normal, 
namely, |ai| > 1 0,2 1 • Among solutions Al— A6 and SI— S5, it is natural to select solutions Al and 
SI since the magnitude of their ai and the ratio 1 02/01 1 are closer to those in solutions A and S 
(see Table IX]) . As stressed before, a large Tp is definitely needed in order to accommodate the 
data of Df — > p + i]. 

B. SU(3) breaking and final-state interactions 

As discussed before, based on the exact flavor SU(3) symmetry and our best solutions to the 
topological amplitudes extracted from CF charm decays, the predicted branching fractions of SCS 
D decays are too large for tt + it~ , 7r°7r° and too small for K + K~ , 7r + r/ and 7t + t/. In the following, 
we shall first examine SU(3) breaking effects in color-allowed and color-suppressed tree amplitudes 
within the factorization approach (for an earlier study of SU(3) violation in charm decays, see 



491]). It turns out that while part of the SU(3) breaking effects can be accounted for by SU(3) 
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symmetry violation manifested in T and C amplitudes, in some cases such as the ratio R = T(D° — > 
K + K~)/T(D° — > tt + tt~), SU(3) breaking alone in spectator amplitudes does not suffice to account 
for R. This calls for the consideration of SU(3) violation in the IF-exchange amplitudes. 



1. SU(3) breaking 

To illustration the effect of SU(3) symmetry violation in external and internal H^-emission 
amplitudes, we consider the following modes: D + — > tt + tt° ,ir + rj, tt + ?/ and D° — > K + K~, tt + tt~ 

D+ -)■ 7T+7T 

Since the isospin of the tt + tt° state is I = 2, this channel does not have nontrivial final-state 
interactions. Indeed, it does not receive weak annihilation contributions. This means that the 
short-distance approach should suffice to describe this decay. The factorizable amplitudes read 

Cf m = ^ aMml - ml)F^(ml). (28) 

Hence, 

ml -ml F D «(ml) U_ F^(mp 

T ml -m\ F DK (mlY C f K F^(m 2 K )' { 

where use of Eq. (f22l) has been made. Numerically, we find T' nir /T = 0.96 and C' n7T /C = 0.78. 
Therefore, 

A{D + tt+tt ) = -7=^^(0.96 T + 0.78 C). (30) 
v2 

Because of the less destructive interference between color-allowed and color-suppressed amplitudes 
due to SU(3) breaking effects, B(D + -»■ tt+tt ) is enhanced from (0.89 ± 0.10) x 10~ 3 to (0.96 ± 
0.04) x 10~ 3 , in better agreement with the experiment. 



The relevant factorizable amplitudes here are 

K Vq = ^aifniml-m^F^iml), 

C\ = ^a 2 f q (m 2 D -ml)F D *(ml), 

CL = %a 2 f s (ml-ml)F^{ml), (31) 



where f q , f s are the decay constants of r\ q and rj s , respectively. The masses of rj q and n s read [23J 



m \ q = -y^(°l m n^75« + m d di^d\r] q ) + -^(0\^GG\r} q ) « ml + -^(0\^-GG\r] q ) 

ml = l(0\m s si l5 s\ Vs ) + j s (0\^GG\7 ]s )^2ml-ml + l(Q\^GG\ Vs ) , (32) 
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where contributions to their masses from the gluonic anomaly have been included. We shall use 
the parameters extracted from a phenomenological fit 21, l^ : <f> = (40.4 ± 0.6)° and 



1 a 

■{0\-^-GG\r] q ) = 0.265 ± 0.010, 



V2f q '4vr 

1.352 ±0.007, f q = U. (33) 



(0\^GG\ Vq ) __f s 



>/2(0\%GG\r, s ) fq 

Assuming that the unknown form factor F^ 9 is the same as , the decay amplitudes are 
modified to 

A(D + -»■ vr+r/) = -=V*LV ud (0.85 T + 0.93 C + 2A) cos <f> - V* s V us 1.28 C sin 0, 
v2 

A{D+ -»• vr+r/) = -7=V£Ki«i(0.85 T + 0.93 C + 2A) sin + V^V US 1.28 C cos 0, (34) 
v2 

where we have assumed A' = A. It follows that B{D+ -> vr+r?) = (2.68 ± 0.29) x 10~ 3 [(3.54 ± 
0.21) x 10~ 3 ] and B(D+ -> vr+r/) = (5.03 ± 0.32) x 10~ 3 [(4.68 ± 0.30) x 10~ 3 ], where the data 
are shown in the squared brackets. Hence the discrepancy between theory and experiment is 
substantially improved. Presumably, a better agreement will be achieved if SU(3) violation in 
the ^-annihilation amplitude is included. For example, if A = 0.7A is taken, we will have 
B(D+ -> vr+? ? ) = (3.02 ± 0.19) x 10~ 3 and B{D + -> vr+r/) = (4.69 ± 0.21) x 10~ 3 , bringing our 
predictions closer to the observed data. 

D° -> ir+Tr-,K + K- 



Experimentally, the rate of D — > K + K~ is larger than that of D° — > tt + tt~ by a factor of 
2.8, while they should be the same in the SU(3) limit. This is a long-standing puzzle since SU(3) 
symmetry is expected to be broken at the level of 30%. Without the inclusion of SU(3) breaking 
effects in the topological amplitudes, we see from Table ITT1 that the predicted rate of K + K~ is even 
smaller than that of 7r + 7r~ due to less phase space available to the former. In the factorization 
approach we have 

T'kk _ f K ml-ml F^jml) 

U ml -ml CK) • 1 > 

Using the form factor parametrization given by Eqs. (|27p and input parameters given in Table ITXl 
and assuming no SU(3) symmetry breaking in E so that E' KK = E' n7T = E, we obtain T' KK /T^ n = 
1.32, implying that B(D° -> K + K~) = (3.4±0.1) x 10~ 3 and B(D° -»• vr+vr-) = (2.1±0.1) x 10~ 3 . 
Therefore, SU(3) symmetry breaking alone in T" is not adequate to resolve the puzzle. This calls 
for the consideration of SU(3) flavor symmetry violation in the VF-exchange amplitudes. We will 
come back to this issue in the next subsection. 

2. Final-state interactions 

We learn from the above few examples that when SU(3) breaking effects in color-allowed and 
color-suppressed amplitudes are included, the discrepancy between theory and experiment for SCS 
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FIG. 1: Contributions to D° — > K°tt° from the color-allowed weak decay D° — > K~ir + followed by 
a resonant-like rescattering (a) and quark exchange (b). While (a) has the same topology as the 
VF-exchange graph, (b) mimics the color-suppressed internal VF-emission graph. 



decays will be significantly reduced. However, in order to get a better agreement, violation of SU(3) 
symmetry in weak annihilation (VF-exchange and VF-annihilation) must be taken into account. 
This is particularly true for the ratio R = T(D° K + K-)/T(D° vr+vr"); SU(3) violation in 
spectator amplitudes accounts for only half of the ratio R. 

Under the factorization hypothesis, the factorizable TU-exchange and VF-annihilation ampli- 
tudes are proportional to 02 and cq, respectively. 8 They are suppressed due to the smallness of 
the form factor at large q 2 = m? D . This corresponds to the so-called helicity suppression. At first 
glance, it appears that the factorizable weak annihilation amplitudes are too small to be consistent 
with experiment at all. However, in the diagrammatic approach here, the topological amplitudes 
C, E, A can receive contributions from the tree amplitude T via final-state rescattering, as illus- 
trated in Fig. Q] for D° -)• K°ir°: Fig. Q^a) has the same topology as the VF-exchange diagram E, 
while H](b) mimics the internal IF-emission amplitude C. Therefore, even if the short-distance W- 
exchange vanishes, a long-distance H^-exchange can be induced via inelastic FSIs. Historically, it 



was first pointed out in 52|] that rescattering effects required by unitarity can produce the reaction 
D° — > K cj), for example, even in the absence of the ^-exchange diagram. Then it was realized 
that this rescattering diagram belongs to the generic IU-exchange topology 0]. 

There are several different forms of FSIs: elastic scattering and inelastic scattering such as 
quark exchange, resonance formation, etc. For charm decays, it is expected that the long-distance 
weak annihilation is dominated by resonant FSIs. That is, the resonance formation of FSI via 
qq resonances is usually the most important one due to the fact that an abundant spectrum of 
resonances is known to exist at energies close to the mass of the charmed meson. Indeed, the 
sizable magnitude of E and its large phase determined from experiment are suggestive of nearby 



In the QCD factorization (QCDF) approach 1501. E and A are proportional to the Wilson coefficients ci 
and C2, respectively, and hence A — (02/01 )E |2oL l5ll]. However, the charm quark is not heavy enough to 
justify the use of QCDF in charm decays. 
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resonance effects. (As noted in [531 ] . weak annihilation in VP systems receives little contributions 
from resonant FSIs. We shall show below that the W^-annihilation amplitude in Df —> tt + lo arises 
from final-state rescattering via quark exchange.) A direct calculation of the resonant FSI diagram 



is subject to many theoretical uncertainties. Nevertheless, as emphasized in 54j, [55J], most of the 



properties of resonances follow from unitarity alone, without regard to the dynamical mechanism 



that produces the resonance. Consequently, as shown in [53j, [54|, the effect of resonance-induced 
FSIs can be described in a model-independent manner in terms of the mass and width of the nearby 
resonances. It was found that weak annihilation amplitudes are modified by resonant FSIs as (see 
'•//•• [3]) 

2iS r i \ f i T\ 2i<5. - - ( & 



E = e + (e** - 1) [e + -j , A = a + (e** - 1) la + - j , (36) 
with 

e 2iS r =1 _i E (37) 

mo — TUR + iT/2 

where the VT-exchange amplitude E and IT-annihilation A before resonant FSIs are denoted by 
e and a, respectively. Therefore, resonance-induced FSIs amount to modifying the VT-exchange 
amplitude and leaving the other quark-diagram amplitudes T and C intact. We thus see that even if 
the short-distance weak annihilation vanishes (i.e. e = and a = 0), as commonly asserted, a long- 
distance W-exchange (VF-annihilation) contribution still can be induced from the tree amplitude 
T (C) via FSI rescattering in resonance formation. 

D° ->■ it + tt~,K + K~ revisited 

We have discussed before that resonant FSIs usually give the most important contributions 
to the weak annihilation topology in the PP modes. For AS = 1 D — > Ktt decays, there is a 
J p = + resonance K% (1950) in the sd quark content with a mass 1945 ± 10 ± 20 MeV and a width 
201 ± 34 ± 79 MeV [3l|. 9 Assuming e = in Eq. ([36]), we obtain 

E = 1.68 x 10~ 6 exp(il43°) GeV, (38) 

which is close to the "experimental" value E = (1.53±o;oi) x 10 ~ 6 exp[i(122±l)°] GeV. Presumably, 
a non- vanishing short-distance e will render the phase of E in agreement with experiment. For 
AS = D — > PP decays, there is a nearby scalar resonance /o(1710) which decays to KK and tttt, 
as depicted in Fig. [2J This long-distance contribution has the same topology as the M^-exchange 
diagram. 

Recent lattice and phenomenological studies indicate that /o(1710) is dominated by the scalar 
glueball component 57, [58]. If /o(1710) is primarily a scalar glueball G, it is naively expected that 



T(G — >• 7T7r)/r(G —7- KK) ~ 0.9 after phase space correction due to the flavor independent coupling 
of G to PP. However, experimentally there is a relatively large suppression of rrn production 



The importance of the i^g(1950) contribution to D° — > K 7r + has been emphasized in 56 1 
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FIG. 2: Long-distance resonant contribution to the color-suppressed tree amplitude of D — > 
KK,tttt through the intermediate state /o(1710), where the blob stands for a transition due to 
weak interactions. This has the same topology as the IV-exchange topological diagram. 



relative to KK in / (1710) decays. The ratio r(/ (1710) -> 7rirl/r(/ (1710) -> KK) is measured 
to be 0.41^o'i7 by BES from J/ip — > >y(K + K~, ir + ir~) decays 59]. To explain the large disparity 
between tttt and KK production in scalar glueball decays, it was first noticed by Carlson et al. 



601 ] , by Cornwall and Soni 6l[ and vitalized recently by Chanowitz 



621 ] that a pure scalar glueball 



cannot decay into a quark-antiquark pair in the chiral limit, i.e., A(G — > qq) oc m q . Since the 
current strange quark mass is an order of magnitude larger than m u and wid, decay to KK is 
largely favored over tttt. However, chiral suppression for the ratio T(G — > tttt)/T(G —> KK) at the 
ladron level should not be so strong as the current quark mass ratio m u /m s . It has been suggested 



631 ] that m q should be interpreted as the scale of chiral symmetry breaking since chiral symmetry 



is broken not only by finite quark masses but is also broken spontaneously. 

In short, if the glueball component of /o(1710) is dominant, then this will explain qualitatively 
why the D° meson decays to K + K~ more copiously then ir + ir~ through the resonant FSIs. 



D° -> K°K° 

It is known that the decay D° — > K°K° proceeds only through the VF-exchange mechanism 
and vanishes under SU(3). Since the short-distance TV-exchange is expected to be small and 
furthermore subject to SU(3) cancellation, this mode can only occur through long-distance W- 
exchange induced by final-state rescattering 64J. The nearby pole contribution from /q(1710) will 
contribute to the decay D° — > K°K°. However, in contrast to D° — > K + K~, this channel is 
prohibited under exact SU(3) symmetry. To see this, Fig. [3] shows some possible weak transition 
between D° and /o(1710). It is clear that for K°K° production, both s and d quarks get involved 
and compensate each other in SU(3) limit due to the cancellation of CKM matrix elements, while 
only the s quark gets involved in K + K~ production. It has been shown in 58f] that in the limit of 
SU(3) symmetry, /q(1710) is composed of a scalar gluonium and a small amount of SU(3) singlet. 
Consequently, the D°— /o(1710) weak transition for K°K° production indeed vanishes when SU(3) 
flavor symmetry is exact. 
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FIG. 3: Possible D° — /q(1710) transition from (a) intermediate states such as K + K ,7r + 7r etc., 
and (b) W-exchange. 



C. Annihilation amplitudes Ap and Ay from final-state rescattering 

As discussed before, we have some difficulties in extracting and understanding the W- 
annihilation amplitudes Ap and Ay in D — > VP decays. The topological amplitude expressions of 
Df —7- -K + p° and Df — > tt + u are given by 

A(Df -> ir + p°) = ±V c * s V ud (A v - Ap), 

A(D+ -> ir+u) = ±=V* s V ud {A v + A P ). (39) 

Naively it is expected that Ap = —Ay. The argument goes as follows. The direct ^-annihilation 
contributions via cs — > W — > ud are not allowed in Df — > tt + uj, p + r], p + rj' decays since the (ud) has 
zero total angular momentum and hence it has the quantum number of 7r + . Therefore, G(ud) = — 
and the final states should have an odd G-parity. Since the G-parity is even for ojh + and odd 
for tt + p°, it follows that the former does not receive any direct W-exchange contribution. Can 
one induce Df — > tt + uj from resonant FSIs? The answer is no because the J = 0, 1 = 1 meson 
resonance made from a quark-antiquark pair uihasanoddGparit, As stressed in 0, the 
even-G state it + lo (also prj and prj') does not couple to any single meson resonances, nor to the 
state produced by the T^-annihilation diagram with no gluons emitted by the initial state before 
annihilation. Indeed, a general consideration of resonant FSIs gives the relations [53I . 

A P + Ay = ap + ay, 

A P - A r v = ap -ay + {e 2 ^ - 1) (a P -a v + ^(C P - C v )^j , (40) 

with the superscript r denoting the annihilation amplitude arising from resonant FSIs. The above 
relation shows that Ap + Ay does not receive any qq' resonance (e.g., the 0~ resonance 7r(1800)) 
contributions. Since the above G-parity argument implies ay = —ap, the decay Df — > tt + uj is 
forbidden, whereas Df — > 7r + p° receives both factorizable and resonance-induced VF-annihilation 
contributions. Experimentally, however, it is the other way around: B(Df — > ir+u) = (2.1 ±0.9) x 
10~ 3 , while B(Df — > 7r + p°) is not seen. 
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FIG. 4: Long-distance contributions to — > tt + uj from the weak decay — > p + r/^ followed by 
quark exchange. Plots (a) and (b) have the same topology as Ay and Ap, respectively. Owing to 
G-parity conservation, /j + r/') cannot rescatter into 7r + p°. 



To resolve the above-mentioned puzzle, we note that there are long-distance final-state rescat- 
tering contributions to Df — > tt + uj allowed by G-parity conservation. A nice example is the 
contribution from the weak decay — > p + rj( ) followed by quark exchange (Fig. The 
rescattering of p + rj^ into 7r + p° is prohibited by the G-parity selection rule. Consequently, 
Ap + Ay = Ap + Ay = 2A P , where the superscript e indicates final-state rescattering via quark 
exchange, and Ap — Ay = A r p — Ay. Since — > p + r] has the largest rate among the CF 
—7- VP decays, 10 it is conceivable that Df — > it + uj can be produced via FSIs at the 10" 3 level 
as its branching fraction. Other processes such as the weak decays Df — > j^W°iif(*)+ followed 



by K^°K^+ —> it + uj have been discussed in [6g, |67|. However, this rescattering process does 



not involve quark exchange. This decay may also proceed through pre-radiation of the u. For 
instance, the meson can dissociate into two meson states such as D^°K^ + and D^ + K^° 
which rescatter strongly to (cs)co while the virtual cs state decays weakly to tt + 

As pointed out in [691 ]. the u meson can be produced from the decay Df — > ir + (f) followed by 
uj — (p mixing. This will be possible if 4> is not a pure ss state and contains a tiny qq component. 
Neglecting isospin violation and the admixture with the p° meson, one can parametrize the u-cp 
mixing in terms of an angle 5 such that the physical oj and 4> are related to the ideally mixed states 
uj 1 = (uu + dd) I \/2 and (f) 1 = ss by 



uj \ ( cos 5 sin 5 
d> \ — sin 5 cos 5 




(41) 



and the mixing angle is about \5\ ~ 3.3° [70[] (see 2l|] for the latest determination of 5). Therefore, 



the production of it + uj through uj — (ft mixing is expected to be 

,Pc(vr0) 



B(D+ -> tt+uj)^ mixing = B(D+ -> ir+cft) ( ^4^) ' tan ' *~ 2 ' 2x 10 " 4 ' ( 42 ) 



As discussed before, from the theoretical point of view, B(Df — > p + if) is at most of order 3%. 
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(a) (b) 

FIG. 5: Manifestation of Fig. [J] as the long-distance i-channel contributions to the IF-annihilation 
amplitude of Df — > tt + uj. 



which is too small by one order of magnitude. Hence, it is not the dominant mechanism for the 
Df — > tt+uj decay. 

Since FSIs are nonperturbative in nature, in principle it is extremely difficult to calculate their 
effects. It is customary to evaluate the final-state rescattering contribution, Fig. HI at the hadron 
level manifested in Fig. [5j One of the diagrams, namely, Fig. [5{b), has been evaluated in 67]. As 
stressed in 711 ] ■ the calculation of the meson-level Feynman diagrams such as Fig. [5] involves many 



theoretical uncertainties. If one naively calculates the diagram, one will obtain an answer which 
does not make sense in the context of perturbation theory since the contributions become so large 
that perturbation is no longer trustworthy. Moreover, because the t-channel exchanged particle is 
not on-shell, a form-factor cutoff must be introduced to the vertex to render the whole calculation 
meaningful. We will leave the calculation of Fig. [5] to a separate publication. 

Finally, we would like to point out that the study of the decays Df — > (p + ,af ,a+)uj will 
help understand the annihilation mechanisms in Df decays. Unlike the tt+uj mode, Df — > p+uj 
can be induced by resonant FSIs and hence it is anticipated to have a larger rate. The recent 
CLEO measurements of Df exclusive decays involving an uj show that the branching fractions of 
D+ -> (tt+tt , vr+vr+vr-, ir + r])uj are (2.78 ±0.70)%, (1.58 ±0.46)% and (0.85 ±0.54)%, respectively 



72]. Assuming the validity of the narrow width approximation, B(a + (980) — > tt + t]) = 1 and 



£>(a^~(1260) — > (37r) + ) = 100% in conjunction with the assumption of B(af — > tt+it+tt ) = B(af — 
7r + 7r 7r ), the branching fractions of Df — > a^"(980)u; and Df — > a^(1260)a; are ready to obtain. 



D. Comments on D ^ VV and D (S, A, T)(P, V) decays 



D^VV 

A handful number of VV modes, such as D + — > K*°p + , D° — > 
K*°p°, K*-p+, K*°K*°, p°p°, cpp , K*°uj and Df p+cp, p + uj, K*+K*°, have been ob- 
served. These decays have a richer structure as they have three partial-wave or helicity states. 
Naively it is expected from the factorization hypothesis that the longitudinal polarization is 
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dominant. However, the Mark III measurement [73|] has indicated that the branching fraction of 
the decay D° — > K*°p° is already saturated by the transverse polarization state, while longitudinal 
and transverse polarizations are comparable in D° — > K*~ p + (see, e.g., [74| for a detailed 
discussion). A Dalitz-plot study of four-body D decays will be helpful to resolve the puzzle. 

D ->■ (S,A,T)(P,V) 

Thanks to the powerful Dalitz plot analyses of multi-body decays of D mesons, the structure of 
many resonances can be probed. Besides the vector meson resonances, many parity-even states such 
as the scalar mesons Oq (980), / (980), / (1370), / (1500), / (1710), a (1450), AT£(1430), the axial- 
vector mesons a x (1260), K 1 (1270), ^(1400) and the tensor mesons / 2 (1270) , a 2 (1320), K £(1430) 
have been studied. Branching fractions of quasi-two-body decays D — > SP, VP, TP with S, A, T 
standing for scalar, axial-vector and tensor mesons, respectively, can be extracted from the Dalitz 
analyses of three-body decays of D mesons. The study of D — > SP opens a new avenue to the 
understanding of the light scalar meson spectroscopy, recalling that the underlying structure of light 
scalar mesons a, k, ao(980) and /o(980) is not well established experimentally and theoretically. 

The topological amplitudes for D — > SP, AP,TP decays have been discussed in 7514771]. There 
are several new features. First, one generally has two sets of distinct external ^-emission and 
internal VF-emission diagrams, depending on whether the emitted particle is a party-even meson 
or a parity-odd one. Let us denote the primed amplitudes T" and C for the case when the emitted 
meson is a parity-even one. Secondly, because of the smallness of the decay constants of parity- 
even mesons except for the 3 Pi axial-vector state, it is expected that \T'\ <C \T\ and \C'\ <C |C|. 
This feature can be tested experimentally. Thirdly, since Kq and the light scalars a, k, f$, ao fall 
into two different SU(3) flavor nonets, one cannot apply SU(3) symmetry to relate the topological 
amplitudes in D + — > foir + to, for example, those in D + — > A'g°7r + . Finally, we notice that the new 
data of quasi-two-body decays have been accumulated in the past years to the extent that they 
warrant a serious theoretical study. 



IV. SUMMARY AND CONCLUSION 

We have studied in this work the two-body hadronic charmed meson decays, including both 
the PP and VP modes. The latest experimental data are first analyzed in the diagrammatic 
approach. The magnitudes and strong phases of the flavor amplitudes are extracted from the 
Cabibbo-favored (CF) decay modes using x 2 minimization. The best-fitted values are then used 
to predict the branching fractions of the singly-Cabibbo-suppressed (SCS) and doubly-Cabibbo- 
suppressed decay modes in the flavor SU(3) symmetry limit. In doing so, we have observed SU(3) 
breaking effects more significant than the two-body hadronic bottom meson decays. 

In Section HH we find from the CF PP modes that 

T = 3.14 ± 0.06, C = (2.61 ± 0.08) e"^ 15 ^ 1 ) ; 
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E = (1.53±g : g|) e ^ 122±2 )° , A = (0.39±g;Si) e^-S) (43) 

for the 77-77' mixing angle <j) = 40.4° inferred by KLOE 2l|. Due to the larger branching fractions 
of K°r] and ir + r]^ modes, the magnitudes of T and C quoted above are larger than previously 
obtained. Also, the T and C amplitudes subtend a strong phase of about 150°. Moreover, the E 
and A amplitudes are roughly perpendicular to each other. 

In the case of VP modes, we find two possible solutions for the Ty, Cp, and Ep amplitudes, as 
given in Table IVTI1 The SCS VP modes D°' + — > 7r 0,+ </> favor the solution with a large Cp (namely, 
solutions A and S in Table IVIII depending on which formula is used to extract the invariant 
amplitude). By assuming that Tp and Ty are relatively real, we then obtain several possible 
solutions for Tp, Cy, and Ey, as listed in Table [VTIIl Again, only one of them (solutions Al and 
SI, with the largest Tp) is in better agreement with the SCS VP modes, such as the Df — > p +r q 
mode. In this case, \(T,C, E)p\ > \(T,C, E)y\. The conclusion does not change much even if 
we make the relative phase between Tp and Ty a free parameter because it turns out to be very 
close to zero. Finally, the Ap and Ay amplitudes cannot be completely determined based on 
currently available data. Instead, we use the non-observation of — > ir + p° to argue that they 
are of similar magnitudes, and use the measured rate of Df — > tt + lo to show that they are much 
smaller in size than the color-suppressed ones. However, this leads to a contradiction with data, as 
B(D+ K°K* + ) > B{D+ -> K*°K+) while the former is dominated by Cy and the latter by C P . 
So a full determination of the Ap and Ay amplitudes still await more precise data on the related 
modes. We conjecture that the currently quoted experimental results for both — > K°K* + and 
Df —?• p + rj' are overestimated and problematic. 

In Section IIIH we compare the sizes of color-allowed and color-suppressed tree amplitudes 
extracted in Section[TI]with the effective parameters oq and 02 defined in the factorization approach. 
The extracted 02, as given in Table [Xl is significantly larger than that obtained from short-distance 
calculations. The ratio (02/01 1 is more or less universal among the D — > Kir, K*ir and Kp modes. 
This feature allows us to discriminate between different solutions of topological amplitudes. The 
relative strong phase extracted from the Kp modes and that from the Kir and K*ir modes are all 
in the vicinity of 180°, as expected from naive factorization. 

Some Cabibbo-suppressed modes exhibit sizable violation of flavor SU(3) symmetry. We find 
that part of the SU(3) breaking effects can be accounted for by SU(3) symmetry violation mani- 
fested in the color-allowed and color-suppressed tree amplitudes. However, in other cases such as 
the ratio R = T(D° — > K + K~)/T(D° — > tt + tt~), SU(3) breaking in spectator amplitudes is not 
sufficient to explain the observed value of R. This calls for the consideration of SU(3) violation in 
the W-exchange amplitudes. Since weak annihilation topologies in D — > PP decays are dominated 
by nearby scalar resonances via final-state rescattering (for example, the resonance ^(1950) helps 
explain the magnitude and phase of E), we argue that the long-distance resonant contribution 
through the intermediate state /o(1710) can naturally explain why D° decays more copiously to 
K + K~ than ir + ir~ through the VF-exchange topology. This has to do with the dominance of the 
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scalar glueball content of /o(1710) and the chiral-suppression effect in the decay of a scalar glueball 
into two pseudoscalar mesons. The same FSI through the /o(1710) pole contribution also explains 
the occurrence of D° — > K°K°. However, in contrast to D° — > K + K~ , this channel vanishes under 
SU(3) symmetry owing to the cancelation of CKM matrix elements V* s V us and V* d V u d and the fact 
that /o(1710) is composed of a scalar gluonium and a small amount of SU(3) singlet in the SU(3) 
limit. 

The CF decays Df — > p + tt°, P°tt + , tt + uj, p + oj proceed only through the ^-annihilation topol- 
ogy. Owing to the G-parity selection rule, Df — > tt + oj does not receive contributions from the 
short-distance VF-annihilation and resonant FSIs. Nevertheless, it can proceed through the weak 
decays Df — > p + i]^ followed by the final-state rescattering of p + rj^" > into tt + uj through quark 
exchange. In contrast, Df — > p + uj can be induced by resonant FSIs and hence it is anticipated to 
have a larger rate. 

Finally we have made some brief remarks on other hadronic decays such as D — > VV and 
D — > (S, A, T)P, the detailed investigation of which is left to a future work. 
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